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Abstract. We used the Atacama Pathfinder Experiment (APEX) 12 m telescope to observe the Jk^k^ - 3o3 ^ 2o2, 322 ^ 
221, and 32i 22o transitions of para-H2C0 at 218 GHz simultaneously to determine kinetic temperatures of the dense gas 
in the central molecular zone (CMZ) of our Galaxy. The map extends over approximately 40'x8' (~ 100x20 pc^) along the 
Galactic plane with a linear resolution of 1.2 pc. The strongest of the three lines, the H2CO (3o3 2o2) transition, is found 
to be widespread, and its emission shows a spatial distribution similar to ammonia. The relative abundance of para-H2C0 is 
0.5-1.2x10"^, which is consistent with results from lower frequency H2CO absorption lines. Derived gas kinetic temperatures 
for individual molecular clouds range from 50 K to values in excess of 100 K. While a systematic trend toward (decreasing) 
kinetic temperature versus (increasing) angular distance from the Galactic center (GC) is not found, the clouds with highest 
temperature (r^n > 100 K) are all located near the nucleus. For the molecular gas outside the dense clouds, the average kinetic 
temperature is 65+10 K. The high temperatures of molecular clouds on large scales in the GC region may be driven by turbulent 
energy dissipation and/or cosmic-rays instead of photons. Such a non-photon-driven thermal state of the molecular gas provides 
an excellent template for the more distant vigorous starbursts found in ultraluminous infrared galaxies (ULIRGs). 
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1. Introduction 

The Galactic center (GC) region is the closest galaxy core. It 
is characterized by a high concentration of molecular gas lo- 
cated in the innermost few hundred parsec of the Milky Way, 
the central molecular zone (CMZ) (Morris & Serabyn 1996), 
and by extreme conditions like high mass densities, large veloc- 
ity dispersions, strong tidal forces, and strong magnetic fields. 
Therefore it is a unique laboratory for studying molecular gas 
in an environment that is quite diflTerent from that of the Milky 
Way's disk. For a general understanding of the physics involved 
in galactic cores, measurements of basic physical parameters, 
such as molecular gas density and gas kinetic temperature, are 
indispensable. 

In local dark clouds, gas temperatures can be constrained 
by observations of the /= 1 transition of CO, both be- 
cause this transition is opaque and easily thermalized and be- 
cause the emission even fills the beam of a single-dish tele- 
scope. At a distance of 8 kpc (Reid 1993), however, the beam 
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filling factor of CO / = 1 ^ emission is unknown, and the 
GC clouds may also be aflfected by self-absorption. Analysis 
of multilevel studies of commonly observed linear molecules 
like CO, HCN, or HCO^ suffers from a coupled sensitivity to 
the kinetic temperature and gas density, making an observed 
line ratio consistent with both a high density at a low temper- 
ature and a low density at a high temperature. The metastable 
inversion lines of the symmetric top molecule ammonia (NH3) 
are frequently used as a galactic standard cloud thermometer 
(Walmsley & Ungerechts 1983; Danby et al. 1988). Radiative 
transitions between K-ladders of NH3 are forbidden, and there- 
fore the relative population^depend on the kinetic temperature 
of the molecular gas rather than its density. However, the frac- 
tional abundance of NH3 varies between 10"^ in hot cores and 
10"^ in dark clouds (e.g. Benson & Myers 1983; Mauersberger 
et al. 1987). Furthermore, NH3 is extremely aflTected by a high 



^ NH3 energy levels have quantum numbers (J,K), where J denotes 
the total angular momentum and K its projection on the molecule's 
axis. For metastable levels, J -K. 
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UV flux and tends to show a characteristic "concave" shape in 
rotation diagrams, either caused by a variety of layers with dif- 
ferent temperatures (e.g., due to shocks) or by the specifics of 
the collision rates. Symmetric top molecules such as CH3C2H 
and CH3CN are not widespread and their emission is very faint 
(Bally et al. 1987; Nummelin et al. 1998). Therefore we should 
look for a widespread symmetric or slightly asymmetric top 
molecule that is more favorable for spectroscopic studies to de- 
rive the kinetic temperature of the entire molecular gas. 

Formaldehyde (H2CO) is such a molecule. It is truly ubiq- 
uitous. Wootten etal. (1978) suggested that fractional H2CO 
abundances decrease with increasing density of the gas in star 
forming regions of the Galactic disk. However, this may be 
due instead to decreasing volume filling factors with increas- 
ing density (e.g., Mundy etal. 1987). Unlike for NH3, varia- 
tions in the fractional abundance of H2CO rarely exceed one 
order of magnitude (Johnstone et al. 2003). To give an exam- 
ple: the H2CO abundance is the same in the hot core and in the 
compact ridge of the Orion nebula, whereas the NH3 hot core 
abundance surpasses that of the ridge by about two orders of 
magnitude (Caselli et al. 1993; Mangum et al. 1993). We also 
note that for the starburst galaxy M 82, rkin(NH3) ~ 60 K, while 
for the bulk of the molecular gas rkin(H2C0) ~ 200 K (WeiB et 
al. 2001; Mauersberger et al. 2003; Miihle et al. 2007). 

The relative populations of the K^, ladders of H2CO (see 
Fig. 1 for an energy level diagram) are almost exclusively de- 
termined by collisional processes (Mangum & Wootten 1993). 
Therefore, line ratios involving diff'erent ladders of one of 
the subspecies, either ortho- or para-H2C0, are good tracers 
of the kinetic temperature (Mangum & Wootten 1993; Miihle 
et al. 2007). The energy levels above the ground (para)-H2C0 
state are 10.5 and 21.0 K for the lower and upper states of 
H2CO (3o3 ^ 2o2), and 57.6 and 68.1 K for H2CO (32i ^ 
220) and (322 '^21) (see Fig. [T] for all H2CO energy lev- 
els under 200 K). Therefore, the line ratios are sensitive to 
gas kinetic temperatures less than 100 K, and the uncertainty 
in gas temperature is relatively small for a measured line ra- 
tio at Tkin < 100 K (also see Mangum & Wootten 1993) in the 
case of optically thin H2CO emission. At higher temperatures, 
the H2CO 218 GHz transitions are less ideal because small 
changes in the ratios yield significant changes in Tkin, so that 
the H2CO J = 5-4 transitions at -364 GHz are then becom- 
ing better tracers (Mangum & Wootten 1993). The Jk^k^ = 
3o3 2o2, 322 221, and 32i 220 transitions of para-H2C0 
stand out by being close in frequency. With rest frequencies 
of 218.222, 218.475, and 218.760 GHz, respectively, aU three 
lines can be measured simultaneously by employing a band- 
width of 1 GHz. In this way, the inter-ladder line ratios H2CO 
322 22i/3o3 2o2 and 32i 22o/3o3 ^ 2o2, are free of un- 
certainties related to pointing accuracy, calibration errors, and 
diff'erent beam widths. In this paper, we therefore present ob- 
servations of the H2CO line triplet at 218 GHz to study the gas 
kinetic temperatures of Galactic center clouds. 

2. Observations and data reduction 

Simultaneous measurements of the Jk^k^ = %3 2o2, 
322 221, and 32i 220 transitions of para-H2C0 (see 



^ and Fig. 1) were obtained with the Atacama Pathfinder 
Experiment (APE^^ Glisten et al. 2006) 12 m telescope 
located on Chajnantor (Chile) between 2010 April and 
2010 June. We used the APEX-1 receiver, operating at 
211-270 GHz, which employs a superconductor-insulator- 
superconductor (SIS) mixer with a typical sideband rejec- 
tion > lOdB. As backends we used a fast fourier transform 
Spectrometer (FFTS, Klein et al. 2006), which consists of two 
units with a bandwidth of 1 GHz each and a channel separation 
of 244 kHz. The full-width-at-half-maximum (FWHM) beam 
size was approximately 30" in the observed frequency range 
and the typical pointing error was ~ 3''. The main beam effi- 
ciency and the forward efficiency were 0.75 and 0.97, respec- 
tively. 

We used the on-the-fly observing mode measuring 4'x4' 
maps in steps of 9'' in both right ascension and declination, 
0.8 second integration time per position, and one OFF position 
measurement after every two map rows, i.e., after about one 
minute of observing time. The surveyed area is 302 square ar- 
cmin and its dimension is roughly 40' x8' along the Galactic 
plane. The total observing time was about 41 hours. 

The data were reduced with the CLASS softwar^B We first 
excluded data with a high noise level due to distorted base- 
lines. The spectra were resampled in steps of 15'' and smoothed 
to a velocity resolution of 2 kms"^ The final maps comprise 
4825 points, corresponding to 4825 spectra for each transition. 
To optimize signal-to-noise ratios (SNRs) of integrated inten- 
sity and channel maps, we had to determine the valid velocity 
ranges for the spectra, especially for those cases where the lines 
were weak. To avoid noise from channels without significant 
emission, we first created a high SNR spectrum by smoothing 
all spectra within a 60"x60" box to obtain a masking spectrum, 
for which 1 and were assigned for the channels with SNRs 
higher and lower than 5, respectively. Different line windows 
were automatically deteremined to cover the emission from 
different positions, and the baseline was removed before cre- 
ating the masking spectrum. Then, the final spectra were cre- 
ated by multiplying the raw spectrum by the masking spectrum. 
Because all lines were observed simultaneously and because 
the H2CO (3o3 2o2) data have the best SNR, the masking 
spectrum from this data set has also been applied to the other 
transitions by assuming the same velocity ranges for the emis- 
sion from these transitions (for the details of this technique, see 
Dame et al. 2001; Dame 2011). 

3. Results 

3.1. General characteristics of the molecular gas 

For the first time, observations of the triple transitions of H2CO 
at 218 GHz have been performed in a large area of the GC. 
Fig [2] (left panel) shows the extended line emission from the 
H2CO (3o3 2o2) transition. Molecular gas, revealed by 

^ This publication is based on data acquired with the Atacama 
Pathnder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 

3 http://www.iram.fr/IRAMFR/GILDAS 
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H2CO (3o3 2o2), shows a similar spatial distribution as am- 
monia (Glisten et al. 1981). All prominent features identified 
in ammonia, e.g., the clouds M-0. 13-0.08 and M-0.02-0.07 
(the molecular clouds labeled with M followed by the galac- 
tic coordinates are identified by Gtisten et al. 1981), are clearly 
detected in H2CO (3o3 2o2), and labeled in the figure. A map 
providing the noise level of the full region is also presented in 
Fig. [2] (right panel). The median noise value for the 4825 spec- 
tra is about 0.1 K (T^) at a velocity resolution of 2 kms"^ In 
the regions of interest where the line emission is strong, the 
noise level is ~ 0.08 K (T^). To show the line emission in more 
detail, we also present integrated intensity maps (Fig. [3|, as 
well as channel maps (see § 33 ), for all transitions in a zoomed 
region in Fig. [2] (left panel). The distribution of the velocity in- 
tegrated line emission of the observed three transitions is quite 
similar except that the weaker two H2CO transitions are not de- 
tected in some of the regions where the 3o3 2o2 transition is 
still observed. 



3.2. Individual spectral lines 



Spectra from some positions of interest are shown in Fig. |B.l 
with overlaid Gaussian fit profiles, and their locations are 
marked in Fig. [2] (left panel). To achieve better SNRs, we aver- 
aged the spectra within a 30''x30'' area to create a new spec- 
trum for each transition. All transitions from a given position 



are presented in the same panel of Fig. |B.l| but with diff'erent 
off'sets along the y-axis. For those positions with clear detec- 
tions of H2CO 322 ^ 221 and 32i 220, fine parameters are 
listed in Table A.l where integrated intensity, J r^b^^v, peak 
main beam brightness temperature, r^b, local standard of rest 
(LSR) velocity, V, and FWHM fine width, Ayi/2, were ob- 
tained from Gaussian fits. In most cases, only one component is 
needed for the Gaussian fits except in the case of H2CO (322 
221). This line is displaced from the CH30H(422 ^ 3i2) tran- 
sition by only 49 km s"\ and we therefore have applied in this 
case two component Gaussian fits. 

In the central nuclear region, line emission is too weak to 
be detected with sufficient signal-to-noise ratio, while H2CO 
(3o3 2o2) emission is detected southwest, ~3 pc away from 
Sgr A* , showing the broad weak line profile belonging to P5 



(-60'^ -45^0 in Fig.B.l 



The line parameters obtained from Gaussian fits do not 
show many peculiarities when one inspects the central veloc- 



ities and line widths in Table A.l as well as the line pro- 



files in Fig. B.l Sometimes, however, there are two velocity 
components, e.g., at offset positions P9 (60'', 120'') and P12 
(150", 225") with respect to the Galactic center, and the num- 
ber of components for their Gaussian fits had to be doubled. 
In particular, for the 44kms"^ component of P12, the H2CO 
(322 221) integrated intensity seems to be significantly larger 
than that of H2CO (32i ^ 220), because the H2CO (322 ^ 221) 
transition is blended by the emission from CH30H(422 ^ 3i2) 
at 218.440 GHz. 



3.3. Molecular line data cube 

In Fig. [4] channel maps with velocity steps of 18 kms"^ 
are presented, allowing us to clearly separate the components 
based on both their positions and velocities. The total inte- 
grated intensity maps are also shown for comparison in the 
panels on the righthand side. Channel maps in steps of 2 km s"^ 
are presented in the Appendix (Fig. B.2|B.4| ) to show more de- 
tail. 

We choose the H2CO (3o3 2o2) data to describe the in- 
dividual molecular concentrations. H2CO 3o3 2o2 emission 
mainly ranges from -27 to 81 kms"^ In the velocity range 
[-27,-9], there are two prominent features, the -15 kms"^ 
cloud M 0.02-0.05 at P12, and the southernmost cloud at PI. 
The latter is part of the 20kms-^ cloud M-0. 13-0.08. There 
is also weak H2CO 3o3 2o2 emission at P5 close to Sgr A*, 
which traces the southwestern lobe of the circumnuclear disk 
(CND). 

Within the velocity range [-9,9], the bulk of the 20kms"^ 
cloud appears in the south and its size is about 7pcxl5pc. 
There are three small cores at P9, Pll, and PI 2 with sizes of 
~ 1 to 2pc in the northeast and one small core at PI 7 even 
farther away from Sgr A*. The peak of the 20kms"^ cloud 
moves from south (-75", -390") relative to the Galactic cen- 
ter to north (-30", -210"), increasing the velocity from -27 to 
-h27 kms"^ 

Within the velocity range [9,27], there is a dense concentra- 
tion at P7 with a size of ~ 2.7 pcx5.4pc. Extended weak emis- 
sion is detected around PI 9 in an irregular morphology with a 
size of ~ 2pcx9pc. The northernmost core, M 0.25-1-0.01, be- 
gins to appear at P22. This concentration is not fully covered 
by our observations. 

In the velocity range [27,63], the line emission shows a 
complex morphology. The prominent features are the 50 km s"^ 
cloud M-0.02-0.07 around P6, P7, and P8, two compact con- 
centrations M 0.07-0.08 at PI 3 and M 0.1 1-0.08 at PI 5, an ex- 
tended region of weak emission associated with M 0.06-0.04 
around P12 and P14, a concentration M 0.10-0.01 around P16, 
and the northernmost core M 0.25 -hO. 01 at P22. In the velocity 
range [27,45], the gas in the southeast close to Sgr A* appears 
to connect the 20 kms"^ and 50 kms"^ clouds. 

In the extreme velocity range [63,81], the 50 kms"^ cloud 
moves to the west by -60" with respect to [27,63], and peaks 
at PIO, the edge of the cloud. In addition to features identified 
in the previous velocity ranges, there are three clumps around 
PI 8, P20, and P21, with sizes of about 3 pcx6pc. 

4. Discussion 

4.1. Formaldehyde column density and abundance 

In the following we derive H2CO column densities and abun- 
dances. Assuming the line emission is optically thin and the 
contribution from the cosmic microwave background is negli- 
gible, the H2CO column density, N(Jk^kJ, in an upper state 
-^KaKc can be obtained by 



3k 



-J 



Tmb dv , 



(1) 



4 



Y. Ao et al.: The thermal state of molecular clouds in the Galactic Center 



(for the equations in this subsection, see Mangum & Shirley 
200^, where h is Planck's constant, k denotes Boltzmann's 
constant, fi is the dipole moment, v the frequency of the transi- 
tion, Tex the excitation temperature, S the line strength, Tmb the 
main beam brightness temperature, and J dy the integrated 
line intensity for the transition Jk^k^ ^ C-^- l)KaKc-i • H2CO is a 
slightly asymmetric top molecule, and its line strength, S , can 
be approximately calculated as for a symmetric top molecule 
by 



7(27 + 1) 



(2) 



The total column density, Motai, is related to the column density, 
N(Jk,k,), in the upper state 7k,k, by 

Motal Z ,E(Jk^ 



N(Jk,kJ gjgKgi kT^x 



(3) 



where gj (=2J+1) is the rotational degeneracy, marks the 
K degeneracy, gi the nuclear spin degeneracy, £(/KaKc) the en- 
ergy of state /KaKc above the ground level, and Z the partition 
function. The partition function Z can be calculated by 



CO / j-K^ + i^Kc<J 

/=0 Ka^O Kc^J-Ka 



gjgKgi exp( — ). 

f^J^ ex 



(4) 



For para-H2C0, gK=l, gi=^, and can only be even. Here we 
include 41 levels with energies above ground state up to 286 K 
and assume the same excitation temperatures for all transitions 
(i.e., local thermodynamical equilibrium, LTE) to estimate the 
partition function. Substituting units and parameters for H2CO 
(3o3 2o2) in Equation (3), and assuming an excitation tem- 
perature of 10 K, the total para-H2C0 column density, Motai, 



IS 



Motal = 1.32x10'^ ^mb(3o3 2o2)<ivcm 



(5) 



^03 



2o2) dy, is in 



where the integrated line intensity, J rmb(3o 
units of K kms"^. 

For the Sgr A* complex with a size of 30 pc in diame- 
ter, corresponding to ~ 13 arcmin, the estimated molecular gas 
mass is -0.4x10^ M© (Kim etal. 2002), yielding an aver- 
age H2 column density of 2.6x 10^^ cm"^. In this region, 
2161 formaldehyde spectra were obtained with a spacing of 
15'', and the average spectrum has an integrated intensity of 
10.5 K km s"^ on a r^b scale. The derived averaged para- 
H2CO column density is 1.4x10^^ cm"^ assuming an excita- 
tion temperature of lOK. Adjusting the excitation temperature 
in the range of 5 to 40 K, the resulting total para-H2C0 column 
density will decrease with increasing Tgx from 5 K to around 
14 K, because with increasing Tex the populations of the / = 
2 and 3 states become more dominant. Beyond Tex = 14 K the 
resulting column densities will increase with excitation temper- 
ature because then also the / > 3 levels will be populated. For 
Tkin ^ 100 K (see § 43 for large velocity gradient (LVG) mod- 
eling), Tex values ^40 K require densities ^10^ cm"^, which 



^ https://safe.nrao.edu/wiki/pub/Main/MolInfo/column-density- 
calculation.pdf 



are unrealistically high on a large spatial scale (e.g.. Glisten 
& Henkel 1983). Therefore higher Tgx values can probably be 
excluded and the corresponding averaged para-H2C0 column 
density can be constrained to (1.3-3.1) x 10^^ cm"^. The re- 
sulting para-H2C0 abundance is (0.5-1.2) x 10"^. This abun- 
dance agrees with the values found by Glisten & Henkel (1983) 
and Zylka et al. (1992), who used K^ = l ortho-formaldehyde 
K-doublet absorption lines to obtain H2CO abundances of 
~10"^^ to 2x10"^ in the Galactic center region. Thus it is rea- 
sonable to adopt two fixed limiting para-H2C0 abundances of 
^ and 10"^^ for the LVG analysis in § 



10 



4.3 



4.2. Line ratios 

In this survey, the three 218 GHz rotational transition lines of 
H2CO are observed simultaneously at the same angular resolu- 
tion, providing good data sets to derive the H2CO line ratios. 
Before quantitatively determining gas kinetic temperatures, we 
first present H2CO 322 ^ 22i/3o3 ^ 2o2 and 32i 22o/3o3 ^ 
2o2 line ratios in Fig. [5] as probes of gas temperature. The 
line ratio maps are derived from both channel maps and to- 
tal integrated intensity maps, and the ratios are calculated by 
integrating channels where the H2CO(3o3 2o2) line emis- 
sion is detected above 5 cr. Since the excitation conditions for 
322 22i/3o3 2o2 and 32i 22o/3o3 ^ 2o2 are very similar, 
these two line ratio maps should be nearly identical, and indeed 
the differences mainly arise from the CH3OH contamination in 

H2CO 322 ~^ 221. 

As seen in Fig. [5j the H2CO 32i ^2o/%3 2o2 line 
ratio, with a median value of 0.23, varies significantly across 
the mapped region, from about 0.15 at the edge of the clouds 
to -0.35 toward the 20 km s"^ GMC and the compact concen- 
tration M 0.11-0.08 at P15. In case of narrow one-component 
features, the 322 ^ 22i/3o3 2o2 ratios follow the same trend 
as the 321 ^ 22o/3o3 ^ 2o2 ratios. For the 20 kms"^ GMC, 
the line ratio is higher in the south than in the north in the ve- 
locity range [-9,9]. Within the velocity range [9,27], the ratio 
becomes higher at P3 in the north of the 20 km s"^ GMC. The 
molecular clouds around PI 8, PI 9, an d P20 have two velocity 
components as presented in Fig. B.l one at ~ 25 kms"^ and 
another at ~ 78 km s"\ and are characterized by low line ratios 
as clearly shown in Fig. [5](bottom). Higher line ratios tend to 
suggest higher gas kinetic temperatures and vice versa because 
the relative populations of the ladders of H2CO are almost 
exclusively determined by collisional processes (Mangum & 
Wootten 1993). To be more quantitative and to relate the line 
ratios to kinetic temperatures, we need to adopt LVG radiative 
transfer modeling, which is done in the following section. 



4.3. Kinetic temperatures of the Galactic center clouds 

To evaluate gas kinetic temperatures, we selected the positions 



with the Gaussian fits listed in Table A.l (see also Fig. B.l 



To investigate the gas excitation from the H2CO line measure- 
ments, we use a one-component LVG radiative transfer model 
with collision rates from Green (1991) and choose a spheri- 
cal cloud geometry with uniform kinetic temperature and den- 
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sity as described in Mangum & Wootten (1993). Dahmen etal. 
(1998) estimated that the velocity gradient ranges from 3 to 
6 kms"^ pc"^ for Galactic center clouds. Here we adopt two 



fixed para-H2C0 abundances of [para-H2C0] : 
(see § ^ 



4.1 



10"^ and 10 



-10 



, and a velocity gradient of (dv/dr) = 5 ( km s ) pc 



The modeled parameter space encompasses gas temperatures, 
Tkin, from 10 to 300 K with a step size of 5 K and H2 number 
densities per cm^, log^H2' from 3.0 to 7.0 with a logarithmic 
step size of 0.1. According to Green (1991), collisional excita- 
tion rates for a given transition are accurate to -20%. 

We first choose the H2CO 322 '^21 and 3o3 2o2 data as 
input parameters. Although the 322 '^21 transition is blended 
with CH3OH at a few positions, its line emission is stronger 
than that of H2CO 32i 220, and the components can be sep- 



arated in all cases. This even holds for position PI 2 (see § 3.2 ), 
where we have determined the central velocities of the two ve- 
locity components from H2CO 3o3 2o2 to then carry out 
a four-component Gaussian fit to both H2CO 3o3 2o2 and 
CH3OH 422 ^ 3i2. Comparing computed line intensities and 
their line ratios with the corresponding observational results, 
we can constrain the kinetic temperature. The gas density is not 
well known because it is highly dependent on the adopted frac- 
tional abundance, velocity gradient, and filling factors. Here we 
choose a filling factor of unity to fit the data, implying that we 
obtain beam averaged quantities. 

In Fig. |6jtop), an example is presented to show how the 
parameters are constrained by the reduced distribution of 
H2CO line measurements in the T^m-n parameter space. As 
can be seen, observed H2CO 322 '^2i/%3 '^02 line ra- 
tios (the approximately horizontal lines in the diagram) are a 
good measure of Tkin, independent of density, as long as the 
transitions are optically thin. In contrast, the density is poorly 
constrained, as can be inferred by the diff'erent resulting den- 
sities for the chosen limiting fractional H2CO abundances (in 
the specified case of Fig.[6j the range covers a factor of 3-4 if 
assuming a filling factor of unity). In Table [T] we present de- 
rived gas kinetic temperatures averaged over 30'' boxes at the 
positions shown in Figs. |2] and B.l It is worth noting that in 
most cases the two adopted para-H2C0 abundances, diff'ering 
by a factor of 10 (§4.1), cause only a slight change of less than 
lOK in kinetic temperature because the 218 GHz H2CO tran- 
sitions remain optically thin in most parameter ranges except 
for gas densities higher than 10"^-^ cm" 
H2CO] = 10"^ For [para-H2C0] = 10 
is well above 10^ and outside the plotted range of densities of 
Fig. |6jtop). This demonstrates that H2CO is a good molecu- 
lar thermometer and can be used to reliably determine kinetic 
temperatures. 

Gas temperatures range from 50 K in the southern part of 
the 20 km s"^ cloud to above 100 K in the 50 km s"^ cloud and 
the molecular core of M 0.07-0.08. While a systematic trend 
of (decreasing) kinetic temperature versus (increasing) angular 
distance from the nucleus is not found, the clouds with highest 
temperature (T^in > 100 K) are all located near the center. To 
estimate the overall gas temperature on a large scale of ~90pc, 
we mask the dense clouds by clipping all the emission above 
3 cr in Tl of the H2CO (3o3 ^ 2o2) fine within a 60''x60'' 
box and create an averaged spectrum, yielding a gas kinetic 



^ in the case of [para- 
this limiting density 



temperature of 65+10 K for more diff'use molecular gas outside 
of the dense cores in the GC region. 

If we instead use the H2CO 32i 22o/3o3 2o2 line ratio 
to constrain gas properties, the kinetic temperature is somewhat 
more sensitive to the gas density so that this line ratio is not 
quite as good as a thermometer (see Fig.[6jbottom)). Therefore, 
we focus exclusively on the H2CO 322 '^2i/%3 2o2 line 
ratio to derive the gas kinetic temperature in this study. 

High gas temperatures were first deduced from the 
metastable transitions of ammonia (Giisten etal. 1981, 1985; 
Hiittemeister etal. 1993). Using the CO 7-6/4-3 line ratio, 
Kim etal. (2002) report, for the Sgr A complex, a gas kinetic 
temperature of 47 K on a linear scale of 30 pc. From mm- and 
submm-line spectroscopy. Oka etal. (2011) deduced tempera- 
tures of at least 63 K in the CND. The 20 km s"^ and 50 km s"^ 
clouds were studied by Giisten etal. (1981, 1985) with species 
like NH3 and CH3CN and by Mauersberger etal. (1986) in 
the (J,K) = (7,7) metastable inversion line of ammonia, yield- 
ing gas temperatures in the range 80-100 K. These results are 
roughly consistent with the temperatures derived by us from 
H2CO. We emphasize, however, that ammonia may be more 
aff'ected than H2CO by a peculiar molecule specific chemistry 
and that the degeneracy between high T^i^ and low ^(H2) or 
vice versa is difficult to overcome for the mm- and submm- 
transitions from linear molecules. 

A high gas temperature will affect the Jeans masses of 
dense cores (e.g., for a molecular cloud, its Jeans' mass is 

= i'-f )-'-' P-'-' = 1-25 (Wl^r'-' Me, 

where n is volume density and T^in is gas kinetic tem- 
perature. Mj -1.2 M0 for Tkin^lOK and ^=10^cm"^, and 
14-39 Mo for Tkin = 50-100 K and ^ = 10^ cm"^), and may af- 
fect the initial mass function (IMF) of star formation, resulting 
in a top-heavy IMF in the Galactic center region (Klessen et 
al. 2007). Indeed, Alexander etal. (2007) suggest a top-heavy 
IMF to explain the observed ring of massive stars orbiting 
about 0.1 pc around the Galactic center. However, Bartko etal. 
(201 1) do not find evidence of a top-heavy IMF at distances be- 
yond 12'' from Sgr A*. Higher densities in the GC clouds with 
respect to those in the spiral arms due to a higher stellar density 
and strong tidal forces may lead to lower Jeans masses and may 
thus counteract the effect of higher temperatures. 

4.4. Heating mechanisms in the GC: turbulent heating 
or cosmic-ray heating 

What heats the dense gas to high temperatures in the GC? The 
four most common mechanisms for heating the gas in molecu- 
lar clouds are (a) photo-electric heating in photon-dominated 
regions (PDRs), (b) X-ray heating (XDRs), (c) cosmic-ray 
heating (CRDRs), and (d) turbulent heating. In star-forming 
regions, gas can be heated by electrons released from nor- 
mal dust grains or poly cyclic aromatic hydrocarbons (PAHs). 
Gas and dust are thermally coupled in very dense regions 
(^H2>10^ cm"^; e.g., Kriigel & Walmsley 1984). However, the 
gas temperatures derived from H2CO are much higher than the 
fairly uniform dust temperatures of ~ 14-20 K in the dense 
clouds of the GC (Pierce-Price etal. 2000; Garcia-Marm et 
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al. 2011; Molinari et al. 2011). Photons can drive such decou- 
pling only at the very surface of irradiated clouds (~ a few per- 
cent of their total molecular gas mass; see Bradford et al. 2003 
and references therein). Only there might they have a chance 
to dissociate complex molecules, such as H2CO, which then 
necessarily probe much deeper and UV-shielded gas regions 
where photo-electric heating of the gas is no longer dominant. 
Therefore, some other process(es) should exist that are effi- 
cient at directly heating the gas from outside. X-ray heating 
(Maloney etal. 1996), cosmic-ray heating (Glisten etal. 1981, 
1985; Papadopoulos 2010), and turbulent heating (Glisten etal. 
1985; Schulz etal. 2001; Pan & Padoan 2009) are such poten- 
tial heating mechanisms. 

A few other heating processes also deserve to be men- 
tioned. Gravitational heating (Goldsmith & Langer 1978; 
Tielens 2005) is important during the collapse phase of molec- 
ular cloud cores. However, gravitational collapse is a tempo- 
rary phase and star-forming activity, with the notable exception 
of SgrB2 (not covered by our survey), is not vigorous in the 
clouds of the CMZ. Heating of the gas by magnetic ion-neutral 
slip (Scalo 1977; Goldsmith & Langer 1978) would require 
further observations of the ionization fraction and the mag- 
netic fields of the molecular clouds (see, e.g., Ferriere 2009; 
Croker et al. 2010), which is outside the scope of this paper. In 
a highly turbulent environment, the effect of a magnetic field 
on the cloud is to decrease the dissipation of kinetic energy, 
i.e., turbulent heating, and to lead to a lower gas kinetic tem- 
perature. Having already mentioned PDRs, we discuss in the 
following XDRs, CRDRs, and the dissipation of turbulence. 

4.4.1. X-ray heating 

In X-ray dominated regions (XDRs), X-rays photoionize atoms 
and molecules, depositing a significant fraction of the primary 
and secondary electron energy in heat (Maloney etal. 1996; 
Hollenbach & Tielens 1999). Unlike UV photons, hard X- 
ray photons are capable of deeply penetrating dense molecu- 
lar clouds and heating large amounts of gas. The X-ray heating 
rate is given by 



Fx = 1.2x 10-^^( 



( 



r)( 



10^ cm"^ erg cm~^ 



) 



1022 



cm~ 



_0 9 -3 -1 

erg cm s 



(6) 



(Maloney et al. 1996), where n is gas density in units of cm~^. 
Fx is X-ray flux density in units of erg cm"^ s"^ , and the col- 
umn density of hydrogen attenuating the X-ray flux. For the 
dense cores studied here, N is of the order of 10^^ cm"^. In the 
region associated with the Galactic center and encompassing a 
diameter of about 20 arcminutes, the total X-ray luminosity in 
the range 2 to 40 KeV (Muno et al. 2004; Koyama et al. 2007; 
Yuasa etal. 2008; Dogiel etal. 2010) is 6x10^^ erg s-\ yield- 
ing an X-ray flux density of 2.0x10"^ ergcm"^ s"^ 

The gas cooling via gas-dust interaction can be expressed 

as 



g-d 



4 X 10-^^ t'IMt^^ - Td) erg cm'^ s" 



(7) 



(Tielens 2005), where Tkin is the gas kinetic temperature and 
dust temperature. For a gas density range of 10^ to lO^cm"^, 
accounting for the velocity gradient, dv/dr, the line cooling is 
approximated well by the expression 



Agas ~ 6 X 10"^^ n^'^ T^^ dv/dr erg cm 



' s-i 



(8) 



(Goldsmith & Langer 1978; Goldsmith 2001; Papadopoulos 
2010), where dv/dr is the velocity gradient adopted to calcu- 
late the line cooling rate in the LVG radiative transfer models. 
The line cooling is through rotational lines of CO, its isotopo- 
logues i^CO and C^^O, and other species (Goldsmith 2001). If 
the water is highly abundant in the warm clouds of the CMZ, 
the line cooling from the water emission will be important, and 
this will lead to lower gas kinetic temperatures in this section. 
However, this will not drastically change our results, because 
water is unlikely to be that dominant. Furthermore, water va- 
por only affects the cooling rate and not the heating process, so 
that our evaluation of relative efficiencies of different heating 
processes is not seriously affected. 

If X-ray heating dominates the heating process, the gas ki- 
netic temperature can be estimated from the thermal equilib- 
rium 



Tx = Ag-d + Ag; 



(9) 



To simply solve the equation above, we set the dust tempera- 
ture to T^=OK, yielding a minimum Tkin value and its simple 
analytic solution as below 



T]dn - ( 



(16 X 10-^ + 720 n"-^ dv/dr Fx (N/lO^^y^y-^ 



12 dv/dr 



4x 10-4^1-^ 
12 dv/dr 



f^^ . (10) 



Choosing a velocity gradient of 5 km s"^ pc"\ and adopting the 
observed X-ray flux density of 2.0x10"^ ergcm"^ s"\ the de- 
rived gas kinetic temperature is only 1 K. It indicates that the 
observed X-ray flux density cannot explain the high gas tem- 
peratures observed by formaldehyde. Peculiar conditions in the 
CMZ, like a potentially high water vapor abundance, enhanc- 
ing cooling, or a high magnetic field inhibiting the dissipation 
of turbulent motion, would not help to diminish the resulting 
discrepancy. 

To simplify this: If the X-ray flux density were about 500 
times higher than the observed value and line cooling domi- 
nates the cooling process, we could derive an approximate gas 
temperature with 



( 



rkin = 45( 



dv/dr 



104-5 cm-3 5 kms-^pc-i 



-1/3 



500 X 2.0 X 10-3 ergcm-2 s 



1023 cm-2 



)-«-^K. 



(11) 



This equation shows that r^in depends weakly on the gas den- 
sity, velocity gradient, colunm density of hydrogen attenuating 
the X-ray emission, and X-ray flux density. The first three pa- 
rameters cannot change significantly. If X-rays really play an 
important role in heating the gas, an X-ray flux density, about 
three orders of magnitude higher than observed, is required. 
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Such an intense X-ray radiation field may exist in some spa- 
tially confined regions. However, it cannot explain the high gas 
temperatures on the large spatial scales of the GC. Recent ob- 
servations (e.g., Eckart et al. 2012; Nowak et al. 2012) confirm 
that the X-ray emission from Sgr A* shows flares, almost daily, 
by factors of a few to ten times over the quiescent emission 
level, and rarely even up to more than 100 times that level on 
time scales from a few minutes to a few hours. However, this 
is still much less than what is required to explain our observed 
kinetic gas temperatures. 

4.4.2. Cosmic-ray heating 

For the UV-shielded and mostly subsonic dense gas cores (e.g., 
dark clouds), cosmic-ray heating is the major heating process, 
and it may also play an important role in heating the gas in 
the GC because the cosmic-ray flux density is enhanced by the 
supernovae in this region. The cosmic-ray heating rate is given 
by 



CR 



3.2x 10-^^n( 



10 



-17 e-1 



) erg cm s 



(12) 



(Goldsmith & Langer 1978), where fcR is the total cosmic-ray 
ionization rate. If the gas heating in the GC is dominated by the 
cosmic-ray heating, we can obtain the gas kinetic temperature 
from the energy balance equation 



FcR - Ag_d -h Ag 



(13) 



As above in § |4.4.1| we set the dust temperature to r^^OK, 
yielding a minimum Tkin value and the following analytic solu- 
tion as below 



T]dn - (- 



(16x 10-^ -h 768^0-5 dv/drj^^y 



,0.5 



12dv/dr 



4 X 10-4^ 
12dv/dr 



1.5 



2/3 



(14) 



Similarly, we can obtain a simple solution as below if line 
cooling dominates the cooling process 



^ kin 



= 6( 



dv/dr 



-1/3 



104-5 cm-3 5 kms"^pc-i 

^CR .1/3 
-1 ^ 



( 



10- 



K. 



(15) 



If cosmic-rays are the only heating source, gas temperatures are 
constrained by three parameters: gas density, velocity gradient, 
and the cosmic-ray ionization rate. It is reasonable to adopt a 
gas density of lO'^-^ cm"^ (still low enough to yield (almost) 
optically thin H2CO 3o3 2o2 emission) and a velocity gra- 
dient of 5 kms"^ pc"^ for the clouds in the GC. The poorly 
determined parameter is the comic-ray ionization rate. Using 
the assumed parameters, the gas kinetic temperatures from 
Equation (14) are 22, 54, 70, and 122 K for cosmic-ray ioniza- 
tion rates of 10"^^ lO'^^^ 2x10-^^^ ^^d 10"^^ s'^ Thus, if cos- 
mic rays play an important role in heating the gas, a cosmic-ray 
ionization rate of at least 1-2x10"^^ s"^ is required to explain 
the observed temperatures in the GC, which is about three or- 
ders of magnitude higher than in the solar neighborhood (e.g.. 



Farquhar et al. 1994). Such an enhanced flux of cosmic-ray 
electrons is inferred in SgrB2 by Yusef-Zadeh etal. (2007), 
and is interpreted as the main molecular-gas heating source in 
this region. The required high cosmic ray flux of 1-2x10"^^ s"^ 
would lead to an Hi density of about 500 cm"^ (Glisten et al. 
1981). Such a large Hi abundance will cause noticeable 21 cm 
signals, which have indeed been seen in Hi absorption surveys 
toward some GC clouds (e.g., Schwarz et al. 1977; Lang et al. 
2010). For comparison, Bradford etal. (2003) estimate that 
a high supernova rate in the nucleus of NGC 253 results in a 
cosmic-ray ionization rate of 1.5-5.3x10"^"^ s"^ This mech- 
anism may also play an important role in regulating the gas 
in ultraluminous infrared galaxies (ULIRGs) where the cos- 
mic ray energy density may be as high as 1000 times that of 
the local Galactic value or even higher (Papadopoulos 2010; 
Papadopoulos et al. 2011). 

4.4.3. Turbulent heating 

The dissipation of turbulent kinetic energy provides a poten- 
tially important heating source in Galactic astrophysical envi- 
ronments, such as interstellar clouds (e.g., Falgarone & Puget 
1995) and the warm ionized medium (e.g., Minter & Balser 
1997), and in extragalactic environments, such as intracluster 
cooling flows (Dennis & Chandran 2005). The observed large 
velocity dispersion in the GC requires energy input to support 
the turbulence because the dynamic timescale is rather short, 
around 10^ years. This implies a high turbulent heating rate. 
Following Pan & Padoan (2009), the turbulent heating rate is 
given by 



Tturb = rifimne, 



(16) 



where n and mu are the number density and the mass of the 
hydrogen atom, ji is the mean molecular weight, yu=2.35 for 
molecular clouds, and e is the average dissipation rate per unit 
mass. The average dissipation rate per unit mass, is given by 



6 = 0.5</L, 



(17) 



where cr^ is the one-dimensional velocity dispersion and L 
the size of the cloud. Replacing the expression for e in equa- 
tion (16) and substituting units, the average turbulent heating 
rate is 



turb 



= 3.3 X 10 



ncrlL ^ erg cm 



(18) 



where the gas density n is in units of cm"^, the one-dimensional 
velocity dispersion cTy is in units of kms"\ and the cloud size 
L is in units of pc. We can relate the one-dimensional velocity 
dispersion and the observed FWHM line widths by the conver- 
sion o-v=yFWHM/2.355 (Pan & Padoan 2009). 

If turbulent heating dominates the heating process, the gas 
kinetic temperature can be estimated from thermal equilibrium 



turb 



= Ag-d + Ag: 



(19) 



To solve the equation above simply, we set the dust tempera- 
ture to rd=OK, yielding a minimum Tkin value and its simple 
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analytic solution as 



5. Conclusions 



^kin 



( 



(16x 10"^ +7920^^-^ dvl&ralL-^f-^ 



U&vl&r 



-4 ^1.5 



4 X 10-4^ 
12dv/dr 



2/3 



(20) 



Choosing a velocity gradient of 5 kms"^pc"\ an observed 
FWHM line width of 20kms"\ and a cloud size of 5pc, the 
derived gas kinetic temperature ranges between 5 1 and 62 K 
for a range in gas density of 10^ to 10^ cm"^. 

For gas densities ^ < 10^ cm"^ and the assumed parameters 
as above, line cooling dominates the cooling process and the 
solution can be simplified to 



Tkin = 62 ( 



104-5 cm-3 



dv/dr 



5pc 5 kms"^pc-i 

^FWHM 



-1/3 



( 



20 kms 



-1 



)K. 



(21) 



This equation underestimates the gas temperature by less than 
10% in comparison with Equation (20), and shows that the 
temperature depends only weakly on the gas density, cloud 
size, and velocity gradient, but strongly depends on line width. 
Considering the line widths in different clouds, we can use 
Equation (20) to calculate the temperatures and present the re- 
sults in the last column of Table [T] To compare the calculated 
temperatures by turbulent heating with the ones derived from 
the H2CO measurements in § |4.3[ we plot both temperatures in 
Fig. [7] In general, the temperatures agree within the uncertain- 
ties, supporting turbulent heating as a good candidate to heat 
the gas to the high temperatures observed in the GC. High gas 
kinetic temperatures were also deduced from NH3 absorption 
lines toward Sgr B2 by Wilson et al. (1982). They suggested 
that the high gas temperatures were caused by turbulence main- 
tained by shearing forces as a consequence of galactic rotation. 

For turbulent and cosmic-ray heating, we cannot distin- 
guish which mechanism dominates the heating of molecular 
clouds in the GC. With a large interferometer such as ALMA, 
one can try to search for molecular clumps with thermal line 
widths, i.e. with line widths that are dominated by thermal mo- 
tion. If such objects can be found, turbulent heating can be ex- 
cluded because the narrow line widths cannot be explained by 
turbulent heating, and the cosmic-ray heating will then be the 
dominant process to heat the gas to high temperatures. Future 
observations of x{e) = (the average ionization fraction) 
can also help distinguish between these two heating mecha- 
nisms, because high cosmic-ray energy densities will boost this 
fraction, unlike turbulence (Papadopoulos 2010 and references 
therein). 

The special thermal state of the GC clouds may be the av- 
erage state of the molecular ISM in ULIRGs, with a direct im- 
pact on their stellar IMF (Papadopoulos et al. 2011). The high 
temperatures of molecular clouds on large scales in the GC 
region may be driven by turbulent energy dissipation and/or 
cosmic-rays instead of photons. Such a non-photon-driven ther- 
mal state of the molecular gas provides an excellent template 
for studying the intial conditions and star formation for the 
galaxy- sized gas in ULIRGs. 



The /KaKe = 3o3 2o2, 322 ^ 221, and 821 220 transitions of 
para-H2C0 were observed simultaneously with the APEX tele- 
scope, covering an area of roughly 40' x8' along the Galactic 
plane with a linear resolution of 1.2 pc, including the Galactic 
center. The main results from these measurements follow. 

(1) The strongest line of the 218 GHz H2CO triplet, H2CO 
(3o3 2o2), is widespread in the mapped region, and its emis- 
sion shows a morphology similar to ammonia (Glisten et al. 
1981). 

(2) The para-H2C0 abundance is found to be 
0.5-1.2x10"^, which is consistent with previous studies 
of formaldehyde absorption lines at cm-wavelengths in the 
Galactic center region. 

(3) Using LVG models, we can constrain gas kinetic tem- 
peratures to be of about 85 K for the Galactic center clouds, 
ranging from 50 to values above 100 K. While a systematic 
trend of (decreasing) kinetic temperature versus (increasing) 
angular distance from the nucleus is not found, the clouds with 
highest temperature (T^in > 100 K) are all located near the cen- 
ter. Molecular gas outside of the dense cores in the Galactic 
center region is characterized by a gas kinetic temperature of 
65+10 K. 

(4) The high temperatures found in the Galactic center re- 
gion may be caused by turbulent heating and/or cosmic-ray 
heating. Turbulent heating can readily heat the gas to the values 
deduced from H2CO. If cosmic-ray heating dominates the heat- 
ing process, a cosmic-ray ionization rate of at least 1-2x10"^"^ 
is required to explain the observed temperatures. The high tem- 
peratures of molecular clouds on large scales in the Galactic 
center region may be driven by turbulent energy dissipation 
and/or cosmic-rays instead of photons. Such a non-photon- 
driven thermal state of the molecular gas make such clouds 
excellent templates for the starbursts found in ultraluminous 
infrared galaxies. 
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Fig. 1. H2CO energy-level diagram up to 200 K. The H2CO 218 GHz transitions observed in this paper are shown in bold. 
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Right Ascension Offset (crcsec) 



Fig. 2. H2CO (3o3 2o2) integrated intensity map (left) and the noise map (right) observed with the APEX in the GC. Left: 
black contour levels for the molecular line emission (on a scale) are -3, 3, 6, 9, 12, 15, 20, 30, 40, 50 cr (Icr = 0.72 K km s"^). 
The black rectangle shows a smaller region where most of the molecular line emission is detected. The molecular cloud labels 
with M followed by the galactic coordinates are adopted from Glisten etal. (1981). The wedge at the side shows the intensity 
range of the line emission on a scale. The beam size of 30'' is shown in the bottom-left corner. Right: Noise map across the 
region. The wedge at the side shows the range of noise on a scale at a velocity resolution of 2 kms"^. The spectra presented 
in Fig. B.l are shown as the plus symbols with numbers. Sgr A* is the reference position and shown as a white cross symbol. The 
dashed line marks the Galactic plane through Sgr A* . 
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Right Ascension Offset (arcsec) 

Fig. 3. Integrated intensity maps for the different transitions observed in the GC. The shown region corresponds to the rectangle 
(black solid line) in Fig. [2] (left). Black contour levels for the molecular line emission (on a scale) are -3, 3, 6, 9, 12, 15, 20, 
30, 40, 50 cr (Icr = 0.72 K km s"^). The wedge at the side shows the intensity range of the line emission on a scale. The beam 
size is shown at the bottom-left corner of the left panel. The dashed line marks the Galactic plane through Sgr A*. 



Table 1. H2CO LVG results 





[para-HsCO] = 10"^ 


[para-HsCO] = 10-^^ 






T'kin 


T'kin 


^ turb 


Position 


(K) 


(K) 


(K) 


PI (-75", -390") 


50+f 


55+1^ 
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P2 (-69", -330") 


75+10 


80^20 


55 


P3 (-30", -210") 


85+10 


105!23 


57 
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59 


P6 (+135", 0") 


90+25 


100+35 


91 


P7 (+195", 75") 
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86 
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^-^-40 


58 
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45 


P12 (+150", +225"/ 44 km s-^ 




55+25 
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P13 (+315". +315") 
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66 
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55+20 
^^-10 




38 


P15 (+420", +435") 


125+2^ 




64 


P16(+165", +570") 


80+35 




61 


P22 (+390", +1020") 






68 



A velocity gradient of 5kms~i pc~i, a filling factor of unity, and the H2CO 322 ^2\Po3 ^02 line ratios averaged over 30'' boxes were 
adopted in the LVG models to derive the solutions presented in Column 2-3. 

^ The temperatures in Column 4 are calculated with Equation (20) with a velocity gradient of 5 km s~^ pc~\ a typical gas density of 10"^^ cm"^, 
and a cloud size, i.e., turbulent scale, of 5 pc. 




Fig. 4. Selected velocity-integrated maps for the different transitions observed in the GC. Black contour levels for the molecular 
line emission (on a scale) are -3, 3, 6, 9, 12, 15, 20, 30, 40, 50, 70 cr {Icr = 0.48 K km s"^) for the first six columns, and -3, 3, 
6, 9, 12, 15, 20, 30, 40, 50 a (Icr = 0.72 K km s"^) for the last column, respectively. The wedges at the sides show the intensity 
scale of the line emission on a scale. 
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Fig. 5. H2CO 322 ^ 22i/3o3 2o2 (top) and 821 22o/3o3 ^ 2o2 (bottom) integrated intensity ratio map. Ratios are calculated 
when the H2CO 3o3 2o2 line emission is detected above 5 cr. The top and bottom rows should be nearly identical, and the 
difference mainly comes from the CH3OH contamination in H2CO 322 ^ 221. The wedges at the sides show the line ratios. The 
beam size of 30'' is shown at the bottom-left corner. Sgr A* is the origin for the offset coordinates and shown as a cross. 
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Fig. 6. Example of LVG modeling for P4. Top: reduced distribution (mainly vertical contours) for a single-component LVG 
model fit to the H2CO brightness temperatures (black contours, = 1,2,4), as well as H2CO 822 '^2i/%3 '^02 line ratios 
(mainly horizonal contours) as a function of and Ty^n- The solid lines represent the line ratios: 0.05, 0.1, 0.15, 0.2, 0.25, 
0.3, 0.35, 0.4. The red dashed lines show the observed line ratio and its lower and upper limits. The para-H2C0 abundances 
per velocity gradient, [para-H2C0]/(dv/dr), for the LVG models are 2x10"^^ pc(kms-^)-i (left) and 2 x 10"^^ pc (kms"^)"^ 
(right), respectively. In the left panel the lines with a given H2CO line ratio move downwards (lower T^i^) at high density because 
the H2CO lines start to become saturated; this causes intensity ratios for a given r^in to get closer to unity. Bottom: reduced 
distribution for the H2CO 32i 22o/3o3 2o2 fine ratios. The kinetic temperature is sensitive to the gas density so this line ratio 
is a less suitable thermometer. 




Fig. 7. Gas temperatures estimated by turbulent heating versus those derived from the H2CO LVG models with [para-H2C0] = 
gQ^j^ jjj^g shows the relationship for the cases where both temperatures are the same. 
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Appendix A: H2CO line parameters 



Table A. 1 : Line parameters 





Position 
(arcsec, arcsec) 






(K) 


Vlsr 
(kms~^) 


Ayi/2 

(kms-i) 


j Tmhdv 
(Kkms-i) 


p 1 


-75 


-390 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
^20 


0.75(0.04) 
0.13(0.03) 
0.19(0.03) 


2.9(0.4) 
8.0(2.6) 
2.5(1.2) 


35.4(1.1) 
39.8(6.4) 
38.3(3.2) 


26.4(0.7) 

5.1(0.7) 

7.3(0.5) 


P2 


-60 


-330 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


1.95(0.05) 
0.72(0.04) 
0.65(0.03) 


6.5(0.2) 
7.1(0.3) 
6.7(0.2) 


18.8(0.4) 
13.1(0.8) 
14.3(0.7) 


36.7(0.7) 

9.5(0.5) 

9.5(0.3) 


P3 


-30 


-210 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


2.64(0.04) 
0.83(0.03) 
0.76(0.02) 


16.4(0.1) 
18.4(0.2) 
17.7(0.2) 


19.4(0.2) 
18.0(0.4) 
17.7(0.4) 


51.2(0.5) 
15.0(0.3) 
13.5(0.2) 


P4 


30 


-60 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


0.83(0.02) 
0.21(0.01) 
0.21(0.01) 


33.2(0.2) 
34.0(0.4) 
33.5(0.3) 


20.1(0.5) 
19.5(1.4) 
16.1(1.0) 


16.7(0.3) 

4.2(0.2) 

3.4(0.1) 


P6 


135 





H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


1.08(0.04) 
0.37(0.03) 
0.41(0.02) 


46.2(0.3) 
44.6(0.6) 
44.2(0.4) 


30.2(0.7) 
26.0(1.3) 
24.4(0.9) 


32.7(0.7) 

9.6(0.4) 

10.2(0.3) 


P7 


195 


75 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


1.23(0.04) 
0.40(0.02) 
0.36(0.02) 


40.2(0.3) 
37.3(0.6) 
39.8(0.4) 


36.8(0.6) 
31.5(1.4) 
34.9(1.0) 


45.0(0.7) 
12.7(0.5) 
12.4(0.3) 


P8 


135 


90 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


2.07(0.07) 
0.85(0.03) 
0.81(0.02) 


47.0(0.3) 
48.2(0.2) 
47.5(0.2) 


27.4(0.6) 
26.0(0.6) 
23.0(0.4) 


56.5(1.1) 
22.2(0.4) 
18.6(0.3) 


P9 


60 


120 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


0.60(0.02) 
0.16(0.02) 
0.17(0.01) 


49.3(0.3) 
49.3(...)^ 
48.5(0.7) 


28.6(0.6) 
38.2(4.6) 
28.9(1.7) 


17.3(0.3) 

5.9(0.8) 

4.8(0.2) 


P9 


60 


120 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


0.37(0.02) 
0.12(0.02) 
0.11(0.01) 


-4.2(0.4) 
-4.2(...)^ 
-4.7(0.9) 


19.8(0.9) 
16.9(2.3) 
20.2(2.0) 


7.3(0.3) 
2.1(0.3) 
2.1(0.2) 


PIO 


90 


135 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


1.07(0.04) 
0.32(0.02) 
0.29(0.02) 


55.5(0.2) 
54.8(0.6) 
55.0(0.4) 


30.0(0.6) 
31.4(1.6) 
30.1(1.1) 


31.9(0.5) 
10.2(0.4) 
8.8(0.3) 


P12 


150 


225 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


0.56(0.03) 
0.20(0.03) 
0.12(0.02) 


44.4(0.3) 
44.4(...)^ 
42.4(1.3) 


21.5(0.9) 
18.3(2.4) 
19.8(2.9) 


12.2(0.4) 

3.8(0.5) 

2.3(0.3) 


P12 


150 


225 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


0.72(0.03) 
0.29(0.03) 
0.25(0.02) 


-10.9(0.2) 
-10.9(...)^ 
-10.6(0.4) 


15.7(0.5) 
13.5(1.2) 
11.4(1.3) 


11.4(0.3) 

3.9(0.3) 

2.9(0.3) 


P13 


315 


315 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


1.41(0.05) 
0.41(0.04) 
0.40(0.03) 


51.3(0.2) 
52.2(0.7) 
51.5(0.5) 


22.3(0.5) 
21.1(1.8) 
21.5(1.2) 


31.5(0.7) 

8.7(0.6) 

8.7(0.4) 


P14 


165 


330 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


0.99(0.04) 
0.20(0.03) 
0.28(0.03) 


46.7(0.2) 
45.9(1.0) 
46.7(0.5) 


13.6(0.5) 
12.5(3.2) 
12.7(1.2) 


13.5(0.4) 

2.5(0.4) 

3.6(0.3) 


P15 


420 


435 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


1.83(0.06) 
0.68(0.04) 
0.65(0.02) 


52.9(0.2) 
52.7(0.3) 
52.5(0.2) 


21.8(0.5) 
19.6(0.8) 
19.5(0.6) 


39.9(0.7) 
13.4(0.5) 
12.8(0.3) 


P16 


165 


570 


H2CO 3o3 
H2CO 322 
H2CO 321 


2o2 
^221 
220 


1.09(0.05) 
0.33(0.04) 
0.36(0.04) 


51.0(0.3) 
50.7(0.7) 
50.3(0.6) 


20.9(0.7) 
18.9(2.1) 
16.5(1.4) 


22.8(0.7) 

6.3(0.5) 

5.9(0.4) 


P22 


390 


1020 


H2CO 3o3 
H2CO 322 


2o2 
^221 


1.15(0.06) 
0.37(0.04) 


35.7(0.3) 
35.4(0.6) 


22.9(0.8) 
16.2(1.7) 


26.1(0.8) 
6.2(0.5) 



Continued on next page 
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Table A. 1 - continued from previous page 



Position 




^LSR 


Ayi/2 




(arcsec, arcsec) 


(K) 


(kms-1) 


(kms-i) 


(Kkms-i) 




H2CO 321 ^ 220 0.36(0.04) 


35.8(0.6) 


18.2(1.5) 


6.7(0.5) 



^ To obtain Gaussian fits of the H2CO 322 '^21 transition at locations with two velocity components, we adopt the same 
central velocities as determined for the H2CO 3o3 2o2 transition. 

Appendix B: Online material 

B.1. H2CO spectral lines 
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Fig. B.l. Spectra of different positions marked by the numbers in Fig. [2] All transitions from a given position are presented in 
the same panel but with different offsets along the y-axis on a scale. The Gaussian fits are indicated with green lines, and the 
central velocities of the components are shown as dashed vertical lines. 
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6.2. Velocity channel maps 
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Fig. B.2. Channel maps of the H2CO 3o3 - 2o2 emission. Black contour levels for the molecular line emission (on a scale) are 
from 0. 16 K (2cr) in s.pdf of 0. 16 K. Velocity channels range from -27 to +81 km s"^ in s.pdf of 2 km s"^ . The central velocities 
of the channel maps are shown at the top of each panel. The wedge at the side shows the intensity scale of the line emission. The 
beam size of 30" is shown in the bottom-left corner. 




Fig. B.3. Same as Fig. B.2 but for the H2CO 822 - 221 emission. 




Fig. B.4. Same as Fig. B.2 but for the H2CO 821 - 220 emission. 



